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ABSTRACT Light control of motility behavior (phototaxis and photophobic responses) in green flagellate algae is mediated by 
sensory rhodopsins homologous to phototaxis receptors and light-driven ion transporters in prokaryotic organisms. In the pho- 
totaxis process, excitation of the algal sensory rhodopsins leads to generation of transmembrane photoreceptor currents. When 
expressed in animal cells, the algal phototaxis receptors function as light-gated cation channels, which has earned them the name 
"channelrhodopsins." Channelrhodopsins have become useful molecular tools for light control of cellular activity. Only four 
channelrhodopsins, identified in Chlamydomonas reinhardtii and Volvox carteri, have been reported so far. By screening light- 
induced currents among algal species, we identified that the phylogenetically distant flagellate Mesostigma viride showed photo- 
electrical responses in vivo with properties suggesting a channelrhodopsin especially promising for optogenetic use. We cloned 
an M. viride channelrhodopsin, MChRl, and studied its channel activity upon heterologous expression. Action spectra in 
HEK293 cells match those of the photocurrents observed in M. viride cells. Comparison of the more divergent MChRl sequence 
to the previously studied phylogenetically clustered homologs and study of several MChRl mutants refine our understanding of 
the sequence determinants of channelrhodopsin function. We found that MChRl has the most red-shifted and pH-independent 
spectral sensitivity so far reported, matches or surpasses known channelrhodopsins' channel kinetics features, and undergoes 
minimal inactivation upon sustained illumination. This combination of properties makes MChRl a promising candidate for 
optogenetic applications. 

IMPORTANCE Channelrhodopsins that function as phototaxis receptors in flagellate algae have recently come into the spotlight as 
genetically encoded single-molecule optical switches for turning on neuronal firing or other cellular processes, a technique called 
"optogenetics." Only one of four currently known channelrhodopsins is widely used in optogenetics, although electrical currents 
recorded in diverse flagellates suggest the existence of a large variety of such proteins. We applied a strategy for the search for 
new channelrhodopsins with desirable characteristics by measuring rhodopsin-mediated photocurrents in microalgae, which 
helped us identify MChRl, a new member of the channelrhodopsin family. MChRl exhibits several sought-after characteristics 
and thus expands the available optogenetic toolbox. The divergence of the MChRl sequence from those of the four known chan- 
nelrhodopsins contributes to our understanding of diversity in the primary structures of this subfamily of sensory rhodopsins. 
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The genomes of microorganisms from all three domains of life 
encode retinylidene proteins of shared sequence homology 
that are called type 1 rhodopsins to distinguish them from the 
animal visual pigment proteins called type 2 rhodopsins (1). Ini- 
tially, microbial rhodopsins were discovered in haloarchaea but 
were subsequently also found in proteobacteria, cyanobacteria, 
green algae, cryptomonads, dinoflagellates, fungi, and representa- 
tives of some other phylogenetic groups. Functionally, microbial 
rhodopsins act as light-driven ion pumps or light sensors that 
initiate diverse signal transduction mechanisms (2). 

The first eukaryotic microbial rhodopsins with a sensory func- 
tion were cloned from green flagellate algae, in which they serve as 
light receptors for photomotility responses (3-6). Four channel- 
rhodopsins have been identified to date, ChRl and ChR2 from 



Chlamydomonas reinhardtii (3, 7-9) and VChRl and VChR2 from 
Volvox carteri (6, 10). They contain a 7-transmembrane-helix 
(7TM) domain characteristic of type 1 rhodopsins ( 1 ) followed by 
a conserved but more varied extended C-terminal region. Photo- 
excitation of algal sensory rhodopsins leads to generation of the 
photoreceptor currents that depolarize the plasma membrane and 
trigger a sensory transduction cascade that eventually leads to a 
change in the flagellar beating pattern (11, 12). Heterologous ex- 
pression in Xenopus oocytes revealed that their 7TM domains ex- 
hibit a unique property of light-gated ion conductance (7, 8), so 
far not observed in any other protein. To highlight this feature, the 
name "channelrhodopsins" has been coined to distinguish algal 
sensory rhodopsins from other microbial rhodopsins (7, 8). When 
transfected into and expressed in animal cells, e.g., defined sub- 
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populations of rodent brain neurons, channelrhodopsins enable 
their targeted light activation in tissue culture and in living organ- 
isms (13-15). This approach, called "optogenetics," offers tempo- 
ral and spatial resolution superior to that of conventional electri- 
cal or chemical stimulation. Channelrhodopsins have greatly 
benefited research into brain circuitry and pathogenic mecha- 
nisms of neurological and other disorders, and they are candidates 
for future gene therapy applications, such as for restoring vision 
lost to retinal-degeneration diseases. Although today channelrho- 
dopsins as optogenetic tools are most widely used in neurons, the 
range of their application is rapidly expanding to include glial, 
muscle, and embryonic stem cells (reviewed in reference 16) and 
potentially many more cellular systems in which rapid control of 
the membrane potential and/or the intracellular ion concentra- 
tion is desirable. 

Several intrinsic properties of the four known channelrho- 
dopsins limit their application as optogenetic tools (reviewed in 
references 17 and 18). The most widely used, ChR2, has maximal 
spectral sensitivity at 470 nm (8, 19), but excitation at longer 
wavelengths is preferable to minimize light scattering and absorp- 
tion by biological tissues. VChRl is a red-shifted channelrhodop- 
sin variant (10), but it has slower current kinetics than ChR2, 
compromising its performance at moderate-to-high stimulation 
frequencies. Another limiting property is that photocurrents gen- 
erated by all channelrhodopsins in response to a pulse of contin- 
uous light decrease to a plateau level, a process called "inactiva- 
tion." In the most commonly used channelrhodopsin, ChR2, this 
decrease from the peak amplitude can be as large as 80% (8, 20, 
21), which correspondingly decreases the light-induced mem- 
brane depolarization, requiring more intense or longer light 
pulses to induce a desired biological action. This limitation is fur- 
ther aggravated by the low unitary conductance of channelrho- 
dopsins, which is less than that of common ion channels, as esti- 
mated by whole-cell current noise analysis (22, 23). Taking into 
account the rapidly growing demand for optogenetic tools, it is 
highly desirable to improve these characteristics of channelrho- 
dopsins and to expand the toolbox of variants available for a vari- 
ety of experimental systems. The challenge here is that the mech- 
anism of ion conductance by channelrhodopsins remains poorly 
understood, in particular because their crystal structures are not 
yet available. 

Molecular engineering of the known channelrhodopsins by 
site-specific mutagenesis and chimera construction have been ap- 
plied to produce variants with altered properties and to gain in- 
sights into the structure-function relationships within the chan- 
nelrhodopsin molecule (see, e.g., references 21 and 23 to 29). An 
alternative approach is to search for natural channelrhodopsin 
variants in algae. To screen algal species for candidates of new 
channelrhodopsins with desirable characteristics, we used the 
photoelectrophysiological population assay for recording 
rhodopsin-mediated photocurrents (30). We found that photo- 
currents in the scaly green flagellate Mesostigma viride showed 
faster kinetics and more red-shifted spectral sensitivity than those 
of the earlier-studied organisms. M. viride is an alga that attracted 
considerable interest from evolutionary biologists, because early 
ultrastructural studies (31) and later molecular phylogenetic anal- 
yses (32) revealed that this unicellular organism is closer to land 
plants (Streptophyta) than to other green algae (Chlorophyta). 
We cloned a channelopsin-encoding nucleotide sequence from 
this alga and studied its channel activity by expression in HEK293 



cells. Here we show that this new channelrhodopsin variant, 
MChRl, refines our understanding of what structural features are 
essential for channelrhodopsin function. Furthemore, it has a 
more red-shifted spectral sensitivity at neutral pH than the previ- 
ously available channelrhodopsin with the most red-shifted spec- 
tral sensitivity, VChRl, and surpasses VChRl by having faster 
current kinetics and smaller inactivation, all of which makes 
MChRl a promising candidate for potential optogenetic applica- 
tions. 

RESULTS 

Rhodopsin-mediated photoelectric currents in M. viride cells. 

Photoreceptor currents in native M. viride cells were measured in 
cell suspensions under unilateral laser flash excitation. The re- 
sponse was similar to those observed previously in other photot- 
actic flagellates with intrachloroplast eyespots (for a review, see 
reference 33). However, it was significantly faster, with a peak time 
of <0.5 ms and a half-rise time of 100 ju,s (Fig. 1A). 

The contribution of the late (delayed) photoreceptor current 
to the recorded signal, discovered earlier in other green algae (34), 
is much smaller in M. viride (Fig. 1 A). The second component was 
also not found in the fluence response dependence (Fig. IB), 
which was well fitted with a single exponential function that satu- 
rated at about the same light intensities as those of previously 
studied early photoreceptor currents in other flagellate algae. 
However, the curve in M. viride, in contrast to those in other 
flagellates (35, 36), was well fitted without introducing an addi- 
tional low-saturation component, which means that the amplified 
current component (33) did not contribute to the signals. One 
possible explanation is that localization of this component in 
M. viride cells is unfavorable for registration in the unilateral 
mode of our suspension method (for a more detailed discussion of 
geometric foundations of the method, see reference 37). Alterna- 
tively, this component could be relatively small, another argument 
for which is the observation that photomotility responses in 
M. viride were not highly sensitive to light (reference 38 and our 
unpublished observations). Analysis of signal transduction in 
M. viride is beyond the scope of this report, in which we focus on 
the early photoreceptor current that directly reflects channelrho- 
dopsin function in native cells. 

The action spectrum for the photoreceptor current in suspen- 
sions of M. viride cells is shown in Fig. 1C. The spectral sensitivity 
was estimated by calculating the mean value of the current re- 
sponse between 50 /lis and 2 ms. The first 50 /us was omitted, 
because during this time, a fast photosynthetic signal developed 
(30). To completely eliminate the contribution of this signal, the 
traces were further corrected by subtracting the responses evoked 
by a 680-nm laser flash. This wavelength is ineffective for genera- 
tion of channelrhodopsin-mediated currents or motility re- 
sponses in flagellates but is absorbed by the photosynthetic pig- 
ments. The spectrum was corrected for photon density, taking 
into account the logarithmic fluence response dependence 
(Fig. IB). The spectral maximum at 531 nm was determined by 
Gaussian fit of the main peak (the solid line in Fig. 1C) . The 90-nm 
half-bandwidth of the spectrum is typical of retinal proteins, in- 
dicating that a single rhodopsin species is responsible for the fast 
light-induced currents in M. viride. Significant light sensitivity of 
the electric response is observed even beyond 620 nm. This long- 
wavelength action spectrum and the fast kinetics of the photoelec- 
tric current made this organism a promising candidate for cloning 
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FIG 1 (A) Photoelectric response in a suspension of native cells of Me- 
sostigma viride. Unilateral excitation with a 6-ns laser pulse at 530 nm. (B) 
Fluence-response dependence of the photoreceptor current. Excitation was 
with a 520 ± 20-nm light from a photoflash. (C) Action spectrum of the 
photoreceptor current in native M. viride cells. Solid line, Gaussian fit of the 
main peak; maximum, 531 nm; half-bandwidth, -93 nm (for details, see the 
text), rel. u., relative units. 



a channelrhodopsin likely to be useful for optogenetic applica- 
tions. 

M. viride opsin sequence. The opsin sequence cloned from 
M. viride (593 residues) consists of a 7TM domain and an ex- 
tended C-terminal domain, as is characteristic of other known 
channelrhodopsins. Most of the residues known to form the reti- 
nal binding pocket in bacteriorhodopsin (BR) are conserved. 
However, the overall homology of this sequence to other channel- 
rhodopsins is lower than between any two of them: its 7TM do- 
main shows 38% identity and 52% similarity with ChRl, 35% 
identity and 51% similarity with ChR2, 36% identity and 51% 
similarity with VChRl, and 37% identity and 53% similarity with 
VChR2. The genomes of C. reinhardtii and V. carteri are com- 
pletely sequenced, and each contains two channelopsin genes. 
Complete genomic information is not available for M. viride, and 
there is a possibility that its genome also contains more than one 
channelopsin. The M. viride protein cannot be definitively distin- 
guished as channelrhodopsin 1 or 2 on the basis of sequence ho- 
mology since it is only slightly closer to the two channelrhodopsin 
1 sequences than to the channelrhodopsin 2 sequences (Fig. 2A). 
Also, this sequence lacks two structural features that differentiate 
ChRl/VChRl from ChR2/VChR2. First, there is no Glu residue 
corresponding to Glu87 (ChRl numbering), which is implicated 
by mutagenesis as required for the pH-dependent spectral shift 
characteristic of ChRl/VChRl (21). Also, the positions of Tyr226 
(ChRl) and Asnl87 (ChR2), conserved in VChRl and VChR2, 
respectively, are two of the molecular determinants of spectral 
sensitivity, desensitization and kinetics, by which ChRl differs 
from ChR2 (24). In the M. viride protein, this site is occupied by a 
Trp residue. 

In the two known channelrhodopsin pairs, the more red- 
shifted one is defined as channelrhodopsin 1. The spectral sensi- 
tivity of M. viride channelrhodopsin is more red-shifted than that 
of any previously known species, including the most red-shifted 
VChRl (see below). Therefore, we named it MChRl. 

Unexpectedly, several residues in the 7TM region previously 
considered defining features of the channelrhodopsin family 
based on the four known homologs are not conserved in MChRl, 
as described below under "Mutagenesis studies." The C-terminal 
domain of MChRl shows even lower homology than the 7TM 
domain to the C-terminal domains of other channelrhodopsins. 

Channel activity in HEK293 cells. The light-gated channel ac- 
tivity of Mesostigma rhodopsin was tested in HEK293 cells (Fig. 3 ) . 
It has been reported that the C domain of channelrhodopsins did 
not influence their function in heterologous systems (7,8), so here 
we studied only the 7TM domain of MChRl (331 residues). As 
with the other known channelrhodopsins, the photocurrent gen- 
erated by MChRl under sustained illumination decayed from a 
peak to a plateau, a phenomenon called inactivation (Fig. 3A) . The 
relatively small inactivation of MChRl -generated currents 
showed that this protein is functionally closer to ChRl than to 
ChR2, as the former exhibits a lower degree of inactivation (21). 
The amplitude and the sign of the photocurrent depended on the 
holding potential, as is typical of other channelrhodopsins, with 
the reversal potential close to 0 (Fig. 3B). The mean peak value of 
MChRl photocurrents at the saturating light intensity and a 
— 60-mV holding potential was 276 ± 97 pA (mean ± standard 
error of the mean [SEM]; n = 12), which is comparable to that 
reported for VChRl in neurons (10) and to our own measure- 
ments of VChRl in HEK cells. The plateau amplitude of MChRl 



May/June 2011 Volume 2 Issue 3 eOOl 15-1 1 



Bio' mbio.asm.org 3 



Govorunova et al. 



A VChR2 



B 



MChRl 

ChRl 

ChR2 

VChRl 

VChR2 



VChRl 



MChRl 



ChR1 



ChR2 




MEWffiffTC 

EttmTC 

ev™tc 

GEVimTC 

gev\ottc 



PVILI 


LS 


J 




E 


PVILI 


LS 


N 


LTGLAN 


PVILI 


LS 


N 


mats 


K 


PVLLI 


LS 


N 


m 


I 


PVILI 


LS 


N 


s 


L 



FIG 2 (A) Phylogenetic tree of the channelrhodopsin family constructed by the neighbor-joining method; (B) alignment of channelrhodopsin sequences 
(regions of predicted helices B and C of the 7TM domain) . The conserved Glu residues are shown in green, the conserved His residue at the proton donor position 
in BR (Asp96) is shown in blue, other conserved residues of the retinal binding pocket are shown in yellow, and important residues that are not conserved in 
MChRl are shown in red. 



photoresponses saturated at lower light intensities than the peak 
value (Fig. 4), as in other channelrhodopsins (20, 39), and the 
degree of inactivation increased with the light intensity. When 
successive pulses were applied with a variable dark interval, the 
amplitude of the peak recovered, with a time constant of -3.5 s at 
pH 7.4, which was faster than that in other channelrhodopsins 
(20, 39, 40). However, the level to which the amplitude recovered 
was only about 90% of that elicited by the first flash, even after 
minutes of dark adaptation. This suggests that in addition to this 
short-term process, another, long-term recovery might take place. 
Interestingly, the same behavior of channelrhodopsin-mediated 
electrical responses was observed earlier in native algal cells (34). 

Estimation of the absorption properties of a receptor respon- 
sible for a photobiological response is usually done by action spec- 
troscopy. For the action spectrum to faithfully reflect the absorp- 
tion spectrum of the pigment, the response needs to be measured 
far from saturation, ideally in the linear range of the intensity 
response curve. Also, the action spectrum should be measured in 
response to a short stimulus to avoid secondary photochemistry, 
which may produce a spectrally dependent mixture of intermedi- 
ates, affecting the action spectrum. We measured the initial slope 
of the photocurrent (assessed from the mean amplitude of the 
signal recorded during the first 50 ms of illumination with low- 
intensity light) to evaluate spectral properties of the ground state 
of MChRl in HEK cells. This parameter showed a close-to-linear 
dependence on light, even at intensities when both the peak and 
plateau levels saturated (Fig. 4B), and hence intensities could eas- 
ily be corrected for equal photon densities. 

The action spectra of MChRl light-gated channel activity in 
HEK cells measured as outlined above at three values of external 
pH are shown in Fig. 5A. To reveal the position of their maxima, 
the main peaks of the spectra were fitted with a Gaussian function 
(as shown in Fig. 1C for the spectrum in native M. viride cells). At 
all three pH values, the spectral maximum was at 528 nm, which 
corresponded within a 3-nm accuracy to the maximum of the 
spectral efficiency of the pigment in native M. viride cells (Fig. 1C). 

The long-wavelength spectral sensitivity found in MChRl is 
one of the properties highly desirable for optogenetic applications. 
The most red-shifted channelrhodopsin variant reported earlier 
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FIG 3 (A) Photocurrents in HEK293 cells expressing the 7TM domain of 
MChRl recorded by the whole-cell patch clamp method at different holding 
potentials changed in 20-mV steps from —100 mV (the bottom trace) to 
+40 mV (the top trace). (B) Current voltage dependence of the peak photo- 
current. Data are mean values ± SEM of results from four successive scans. 
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FIG 4 (A) Photocurrents in HEK293 cells expressing the 7TM domain of 
MChRl recorded by the whole-cell patch clamp method. Excitation was at 
530 nm for 2 s, with a 30-s dark interval; relative light intensities, starting from 
the bottom trace, were 100%, 20%, 7%, 3%, and 0.5%. (B) Dependence of the 
amplitudes of the peak and plateau currents (filled squares and filled circles, 
respectively) and the initial slope of photocurrent (open squares) on the light 
intensity. The initial slope of the photocurrent was normalized to the peak 
value at the maximal intensity. 

was VChRl from V. carteri (10). We measured the action spectra 
of the light-gated channel activity of VChRl under the same ex- 
perimental conditions as we used for MChRl and found, using the 
same fitting algorithm, that the VChRl spectrum peaked at 
520 nm at neutral pH (Fig. 6B, solid line), which was ~8 nm less 
than that measured in MChRl (Fig. 6A). The 520-nm peak of the 
VChRl spectrum corresponded to the absorption maximum of 
the isolated pigment measured at this pH (6). We found that acid- 
ification of the external medium led to the appearance of a red- 
shifted shoulder in the VChRl action spectrum (Fig. 6B, dashed 
line). The maximal absorption of purified VChRl was reported to 
shift to 540 nm at low pH, indicating conversion of the pigment to 
a protonated form (6). Our results suggest that only partial con- 
version occurred in HEK cells or that the conductance of the pro- 
tonated state was much smaller than that of the 520-nm form. 



0.6- 



6 0.4- 



O 
£Z 
CD 

o 0.0- 



E 

1 0.8- 

CD 

o 



0.4- 



0.2- 



0.0- 



A Jt 




i i i | i i i i | i i 

1 / 

f 

f 
f 


i i | i i i i | i i 

V 

V? 2 

\ 



450 



500 



550 



600 



Wavelength (nm) 

FIG 5 Low-light- intensity action spectra of photoelectric responses in 
HEK293 cells. (A) MChRl at the external pHs of 9.0 (filled squares), 7.4 (open 
circles), and 5.3 (open triangles). As these data were very close, a single B-spline 
solid line was drawn through the average values at each wavelength. (B) 
VChRl at external pHs of 7.4 (filled squares) and 5.3 (open circles). Data for 
each spectrum are mean values ± SEM of results from 6 to 10 successive scans 
in opposite directions obtained on 2 to 3 cells. 



The main difference between the two red-shifted channelrho- 
dopsins was that the current decay after switching off the light was 
much faster in MChRl than in VChRl. In both cases, the decay 
kinetics was better fitted with two exponentials. At neutral pH, the 
time constants of both components were smaller in MChRl than 
in VChRl (Fig. 6A). In addition, the relative amplitude of the slow 
component was smaller in MChRl, further increasing the rate of 
the overall decay. Due to the small contribution of the slow com- 
ponent to the decay kinetics in MChRl, here we analyzed only its 
fast component. Its rate strongly increased upon acidification of 
the external medium from pH 9.0 to pH 5.3 (Fig. 6B, filled sym- 
bols), whereas in VChRl it only slightly increased (Fig. 6B, open 
symbols). However, the rate of the slow component in VChRl 
showed a stronger dependence on pH than the fast component 
(data not shown). 

In view of possible optogenetic applications, we compared the 
responses of MChRl and VChRl to high-frequency stimulation. 
At a 25-Hz stimulus frequency, the photocurrent generated by 
MChRl was modulated to at least 50% amplitude, whereas that 
generated by VChRl showed practically no modulation (Fig. 7A). 
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The much higher degree of amplitude modulation in MChRl than 
in VChRl derives not only from the faster current decay but also 
from significantly less inactivation of the photocurrent by a series 
of light pulses (Fig. 7A). Consequently, at the level of 50% ampli- 
tude modulation, MChRl demonstrated an ~15-fold-better fre- 
quency response than VChRl (Fig. 7B). This probably reflects the 
adaptation of unicellular M. viride for a rotation frequency higher 
than that of colonial V. carteri. 

Mutagenesis studies. As shown above, wild-type MChRl ex- 
hibited typical light-gated channel activity upon heterologous ex- 
pression, although not all structural features found in other chan- 
nelrhodopsins are conserved in its primary sequence. To test for 
the functional importance of these features, we generated and 
studied three groups of MChRl point mutants. 

(i) V102E and A116E mutants. All other channelrhodopsins 
contain an array of five Glu residues in the predicted second helix, 
but only three of these are conserved in MChRl (Fig. 2B). It might 



be expected that introducing the missing glutamates (correspond- 
ing to Glu83 and Glu97 in the ChR2 sequence) in the MChRl 
sequence would enhance its channel activity. However, both mu- 
tations resulted in a dramatic decrease in the whole-cell current 
amplitude; the plateau current at the saturating light intensity was 
28.3 ± 8.0 pA (mean ± SEM; n = 6) in the Al 16E mutant and 0.3 
± 0.3 pA (mean ± SEM; n = 8) in the V102E mutant. The inten- 
sity of fluorescence of the enhanced yellow fluorescent protein 
(EYFP) in cells that expressed the mutant sequences did not 
change compared to that of the wild type. Inhibition of the chan- 
nel function of the MChRl molecule was also revealed by a signif- 
icant slowing down of the current kinetics (see Fig. 8 for the Al 16E 
mutant; data for the V102E mutant are not shown). In the V102E 
A116E double mutant (glutamate was introduced at both posi- 
tions), the photocurrents were below the detection limit of our 
measurements. 

(ii) A153/H/R mutant. All known channelrhodopsins so far 
have a His residue at the position of the proton donor in BR 
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FIG 8 Photocurrents generated by the Al 16E and C147A mutants and wild- 
type MChRl (WT) in HEK293 cells. Traces were normalized at the plateau 
level to reveal the differences in the decay kinetics. The increased level of the noise 
reflects a decrease in the absolute current amplitudes. Excitation was at 530 nm for 
0.5 s. Time zero corresponds to the end of the light pulse. 



(Asp96), and its replacement with Arg resulted in an increase in 
the stationary current amplitude in ChR2 (15). When Ala found in 
MChRl in this position was replaced with His or Arg, the current 
was completely suppressed without decreasing the EYFP fluores- 
cence in the membrane. 

(iii) C147A and D175A mutants. It has been suggested that 
residues Cysl28 and Aspl56 in ChR2 form a hydrogen bond be- 
tween helices C and D that is important for channel gating (41, 
42) . Its disruption leads to a dramatic increase in the lifetime of the 
channel's open state (25, 43). The corresponding residues (C147 
and D175) are conserved in MChRl. Replacing either of these 
residues with Ala led to a large decrease in the EYFP fluorescence, 
indicating a decrease in protein expression levels. The current am- 
plitude was also dramatically decreased; the plateau amplitude 
was only 0.7 ± 0.5 pA (mean ± SEM; n = 9) in the C147A mutant 
and 7.3 ± 7.3 pA (mean ± SEM; n = 3) in the D175A mutant. In 
agreement with the ChR2 result, analysis of the current kinetics 
showed similar significant increases in the decay time (see Fig. 8 
for the D175A mutant; data for the C147A mutant are not shown). 

DISCUSSION 

Our results show that screening flagellate species by measuring 
rhodopsin-mediated photoelectric currents in vivo is an efficient 
strategy for selecting potential source organisms for cloning new 
channelrhodopsins with characteristics desirable for optogenetic 
applications. Such screening cannot be performed by measuring 
phototaxis itself, i.e., the behavioral response. First of all, if an alga 
exhibits phototaxis, it is not necessarily mediated by a channel- 
rhodopsin; a well-known example is Euglena, which uses a flavin- 
binding adenylyl cyclase as a receptor for light control of motility 
(44). Second, even in those species in which channelrhodopsins 
do act as receptors for phototaxis, a change in the flagellar beating 
pattern occurs as the final result of a not yet fully understood 
amplification cascade (for reviews, see references 45 to 47). Con- 
sequently, the strength of phototaxis in a particular species might 



reflect the efficiency of amplification, rather than the properties of 
the cell's channelrhodopsins. The latter can be probed in algal cells 
only by measuring photoelectric responses. 

MChRl is a channelrhodopsin that we cloned from the flagel- 
late alga M. viride, identified as a candidate species by the above 
approach. The MChRl sequence significantly expands channel- 
rhodopsin diversity, as it deviates from other known channelrho- 
dopsins in several important aspects, thereby refining the se- 
quence criteria for light-gated channel activity. Our results 
showed that two out of five glutamate residues in the putative 
helix B that are conserved in the other four channelrhodopsins 
and that have been suggested to contribute to forming a water- 
containing channel (48, 49) are not essential for the channel func- 
tion, at least in some channelrhodopsins. Similarly, although the 
position of the proton donor (Asp96 in BR) is confirmed to be 
sensitive, the presence of His at this site cannot be considered an 
essential feature of the channelrhodopsin family. On the other 
hand, formation of the predicted hydrogen bond between puta- 
tive helices C and D (Cysl28 and Aspl56 according to ChR2 num- 
bering) is required for proper gating even in such a distant mem- 
ber of the channelrhodopsin family as MChRl. Also, the role of 
Glu87 in pH-dependent color tuning in channelrhodopsins (21) 
is indirectly confirmed by the lack of such tuning in MChRl, in 
which no Glu residue is found at this position. 

Analysis of MChRl properties upon expression in HEK cells 
showed that (i) it exhibits the most red-shifted absorption among 
all studied channelrhodopsins; (ii) the spectral sensitivity of its 
channel activity is pH independent; (iii) its current kinetics is 
significantly faster than that of the other known red-shifted chan- 
nelrhodopsin, VChRl; and (iv) it shows lower inactivation than 
VChRl when stimulated with light of the same wavelength. These 
four qualities are highly desirable for a prospective optogenetic 
tool, especially when high-frequency stimulation is required. 

MATERIALS AND METHODS 

Strain and growth conditions. Mesostigma viride strain CCMP2046 (also 
known as NIES 296) was obtained from the Provasoli-Guillard National 
Center for the Culture of Marine Phytoplankton. Cells were grown at 
25°C in modified Jaworski medium [the concentration of Ca(N0 3 ) 2 was 
increased to 0.17 mM; 0.5 mM KN0 3 , 0.15 juM ZnS0 4 , and 0.04 fiM 
CoCl 2 were added] under a 12-h light (-3,000 lx)-12-h dark cycle. 

Measurements of rhodopsin-mediated photocurrents in M. viride 
cells. Currents were measured with the population assay described in 
reference 30. Two platinum wires immersed in a cell suspension pick up 
an electrical current generated in response to a unilateral excitation flash 
from a Vibrant HE 355II tunable laser (Opotek Inc., Carlsbad, CA), set at 
desired wavelengths. The signal was amplified by a low-noise current am- 
plifier (model 428; Keithley Instruments, Cleveland, OH) and digitized by 
aDigidata 1322A digitizer, supported by pCLAMP 10 software (both from 
Molecular Devices, Union City, CA). 

Cloning and expression of M. viride channelopsin. Total RNA was 
extracted from 500 ml of 1 -week-old culture of M. viride using Trizol 
reagent (Invitrogen, Carlsbad, CA). Synthesis of 3' and 5' rapid amplifi- 
cation of cDNA ends (RACE)-ready first-strand cDNAs and 3' and 5' 
RACE PCR were carried out using the SMARTer RACE cDNA amplifica- 
tion kit (Clontech Laboratories, Takara Bio Company, Mountain View, 
CA) using primers designed according to the opsin sequence fragment 
found in the Taxonomically Broad EST Database (http://amoebidia.bcm 
.umontreal.ca/pepdb/searches/welcome.php). The overlapping RACE 
fragments were combined by fusion PCR, cloned into the pCR2.1-TOPO 
vector (Invitrogen), and fully sequenced. The 7TM domain (encoding 
residues 1 to 331) was inserted between BamHI and NotI sites to replace 
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the VChRl sequence in the pcDNA3. 1/VChRl -EYFP mammalian expres- 
sion vector provided by K. Deisseroth (Stanford University). This vector 
was also used for expression of VChRl in HEK cells. The vector map and 
sequence are available at the EveryVector website. The presence of a flu- 
orescent tag is not expected to affect channelrhodopsin properties, as has 
been shown by quantitative comparison of photocurrents generated by 
YFP-, mCherry-, and myc-tagged ChR2 (50). Point mutations were intro- 
duced using the QuikChange XL site-directed mutagenesis kit (Strat- 
agene, La Jolla, CA). HEK293 cells were transfected using the TransPass 
COS/293 transfection reagent (New England Biolabs, Ipswich, MA). All- 
rrans-retinal (stock solution in ethanol) was added to a final concentra- 
tion of 2.5 p,M. 

Whole-cell patch clamp recording. Measurements were performed 
24 to 48 h after transfection with an Axopatch 200B amplifier (Molecular 
Devices, Union City, CA). The signals were digitized with a Digidata 
1440A digitizer using pCLAMP 10 software (both from Molecular De- 
vices). Patch pipettes with resistances of 2 to 5 Mil were fabricated from 
borosilicate glass and filled with the following solution (in mM): KC1 
(126), MgS0 4 (2), CaCl 2 (0.5), EGTA (5), and HEPES (25, pH 7.2). The 
bath solution contained (inmM) NaCl (150), CaCL, (1.8), KC1 (4), MgCl 2 
(1), glucose (5), and HEPES (10, pH 7.4), unless otherwise indicated. For 
experiments at pH 9, Tris was used in the bath solution instead of HEPES. 
Unless otherwise indicated, the holding potential was —60 mV. Light 
excitation was provided by a polychrome IV light source (TILL Photonics 
GmbH, Grafelfing, Germany) pulsed with a mechanical shutter (Uniblitz 
model LS6; Vincent Associates, Rochester, NY). The light intensity was 
attenuated with the built-in polychrome system or with neutral density 
filters. 

Nucleotide sequence accession number. The nucleotide sequence of 
the opsin sequence cloned from M. viride was deposited in GenBank un- 
der accession no. JF922293. 
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